Paraelectricity in Magnetized Massless QED 
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We show that the chiral-symmetry-broken phase of massless QED in the presence of a magnetic 
field exhibits strong paraelectricity. A large anisotropic electric susceptibility develops in the in- 
frared region, where most of the fermions are confined to their lowest Landau level, and dynamical 
mass and anomalous magnetic moment are generated via the magnetic catalysis mechanism. The 
nonperturbative nature of this effect is reflected in the dependence of the electric susceptibility on 
the fine-structure constant. The strong paraelectricity is linked to the electric dipole moments of the 
particle-antiparticle pairs that form the chiral condensate. The significant electric susceptibility can 
be used as a probe to detect the realization of the magnetic catalysis of chiral symmetry breaking 
in physical systems. 
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Effects of strong magnetic fields in QED have been 
an active research area for many years At present, 
such studies have been reactivated by the observation of 
very strong fields, in the range of 10 12 — 10 16 G, in the 
surface of stellar compact objects. Also, by both the- 
oretical and experimental indications that the colliding 
heavy ions at RHIC and other collaborations can gener- 
ate very strong magnetic fields, estimated to be of order 
eH ~ 2m^(~ 10 18 G) for the top collision, ^s NN 200 
GeV, in noncentral Au-Au impacts at RHIC, or even 
larger, eH ~ 15m^(~ 10 19 G), at future LHC experi- 
ments Q. 

On the other hand, the study of theories of mass- 
less relativistic fermions has recently gained new interest 
in the context of quasiplanar systems, such as pyrolitic 
graphites (HOPG) @, and graphene [f| , because their 
low-energy excitation quasiparticlc spectrum have a lin- 
ear dispersion. The dynamics of those charge carriers 
is described by a "relativistic" quantum field theory of 
massless fermions in 2+1 dimensions 0, @j- 

Massless QED in the presence of a magnetic field ex- 
hibits a peculiar phenomenology. Because of the Landau 
quantization of the fermion's transverse momentum in a 
magnetic field, the dynamics of the lowest Landau level 
(LLL) particles is 1+1-dimensional. This dimensional re- 
duction favors the formation of a chiral condensate, even 
at the weakest coupling, because there is no energy gap 
between the infrared fermions in the LLL and the an- 
tiparticles in the Dirac sea. This phenomenon is known 
as the magnetic catalysis of chiral symmetry breaking 
(MCxSB). The MCxSB modifies the vacuum properties 
and induces dynamical parameters that depend on the 
applied field. This effect has been actively investigated 
for the last two decades ^ n the original studies of 

the MCxSB the catalyzed chiral condensate was 

assumed to give rise only to a dynamical fermion mass. 
Recently, however, it has become clear that besides 
the dynamically generated mass, the MCxSB inevitably 
produces also a dynamical anomalous magnetic moment 
(AMM), because this second parameter does not break 
any symmetry that has not already been broken by the 



chiral condensate and the magnetic field. The dynam- 
ical AMM leads, in turn, to a nonperturbative Landc 
g-factor and Bohr magneton proportional to the inverse 
of the dynamical mass. The induction of the AMM leads 
to a nonperturbative Zeeman effect [T^] ■ 

An important aspect of the MCxSB is its universal 
character. It will occur in any relativistic theory of inter- 
active massless fermions in a magnetic field. The MCxSB 
has been proposed as the mechanism explaining various 
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effects in quasiplanar condensed matter systems 

A drawback of the MCxSB phenomenon is that the 
dynamical parameters (mass and AMM) are extremely 
small even at relatively high fields. Aside from the fact 
that it may be experimentally challenging measuring the 
magnetically catalyzed parameters, there may be cases 
where other competing mechanisms are proposed to ex- 
plain a given magnetic field effect. Consequently, it 
would be nice to have an independent way to experi- 
mentally distinguish the MCxSB from other possibilities. 
One of the main purposes of this letter is to argue that by 
using a weak electric field as a probe, one could obtain, by 
measuring the induced electric polarization, compelling 
evidence in favor or against the existence of MCxSB. 

The electric polarization is found as minus the deriva- 
tive of the electromagnetic free energy with respect to 
the applied electric field. For a weak electric field E, the 
free-energy density can be expanded in powers of E as 



$ = $o - T)E - \E 



(1) 



In a magnetized medium, the coefficients <&o, rj, x, etc., 
may in principle depend on the magnetic field. The sus- 
ceptibility coefficient 77 is different from zero for ferro- 
electric materials [12]. In magnetized QED it is zero, be- 
cause the second term in the r.h.s. of (JTJ) violates parity, 
a symmetry that is not broken neither in massive QED 
nor in the chirally broken phase of massless QED. The 
coefficient x characterizes the lowest order of the system 
dielectric response. It accounts for the electric polariza- 
tion P = xE induced by the applied electric field. This 
term does not break any additional symmetry in (fj}. 
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The electric susceptibility \ in massive QED at strong 
magnetic field can be obtained from the one-loop photon 
polarization operator in an external field configuration 
defined by the gauge = (—Ex3,0,—Hxi,0), which 
implies the existence of uniform parallel magnetic and 
electric fields in the ^-direction. Since the electric field 
is only playing the role of a probe, while we are looking for 
a strong magnetic field effect, wc take the approximation, 
\eE\ < m 2 < \eH\, where m is the electron mass. Hence, 
in the calculation of the photon polarization operator the 
fermion Green's functions depend only on the magnetic 
field, while the effect of the weak electric field can be 
neglected. Thus, the general covariant structure of the 
polarization operator under such conditions is fl3l ] 



IW<?) = 5Z Ka (1(a) \2 ' 
a— 1 ^ ' 



(2) 



where b*jf* are the orthogonal vectors, transverse to q^, 
= -q 2 F tll ,F vp q p + q 2 ± q p , 
b^ = \e^ P xh P qx, bf = % v q v , (3) 



— F pv j H denoting the normalized electromag- 



with i 7 ^ 

netic strength tensor. In ©, n a are scalar coefficients 
depending on the magnetic field. At strong magnetic 
field, where the electrons will all lie in the LLL, only the 
coefficient k 2 is different from zero. Thus, from ^ and 
([3]) we see that only the longitudinal components of Yl pu , 
which in turn depend only on the longitudinal momenta, 
contribute in the strong-field limit. In the static limit, 
<7o = 0, and at small spatial momenta (for a constant and 
uniform electric field, the contribution to the free energy 
of powers of momentum higher than quadratic is zero, 
so we do not need to consider them), the polarization- 
operator coefficient k 2 behaves as [l3|- 14 



K 2 (qo = 0, |"<^| 



0) 



-a\eH\ 

37T771 2 
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(4) 



In this limit, the only component different from zero is 
IIoo = ^2- Its contribution to the electromagnetic free- 
energy density is given by 



$ - *o ~ i J A Q (x 3 )U 00 (x 3 - x' 3 )A Q (x' 3 )dx 3 dx 3 



-XQedE 2 



(5) 



with xqed = a\eH\/3irm 2 , and V the system volume. 
Clearly, high field values \eH\ > m 2 ~ 10 13 G are re- 
quired for the electric susceptibility (xqed) to be sig- 
nificant. Therefore, this polarization effect can only be 
relevant for the astrophysics of compact stars and for 
heavy ion collisions, where such large field strengths can 
exist. It is worthy noticing, on the other hand, that the 
result ([?]) is not valid for m = 0. One can check that if the 
electron mass is taken to zero and the calculations are re- 
peated in the strong-field limit, the electric susceptibility 



becomes zero. In this case, we know that the situation 
cannot be changed by any higher order perturbative con- 
tribution, since the chiral symmetry of the massless the- 
ory is protected against perturbative corrections. Chiral 
symmetry, however, can be broken nonperturbatively via 
MC X SB. 

Our goal now is to find the electric susceptibility in 
the chirally broken phase of massless QED in the pres- 
ence of uniform electric and magnetic fields along the 
^-direction. As already discussed, along with the in- 
duced dynamical mass, the chiral condensate necessar- 
ily produces a dynamical AMM [lfj. Here again, we 
are looking for a strong magnetic field effect and treat 
the electric field as a weak probe. Hence, we assume 
\eE\ < (E ) 2 < \eH\, with E° the dynamical LLL rest 
energy. Accordingly, we neglect E in the fermion propa- 
gator. The one-loop photon polarization operator is 



H-iMu{x,y) = -A-KiaTr[y li G(x,y)~f v G(y,x)] 
where the electron full propagator is given by 



G(x,x') 



(6) 



(7) 



with = J2J ^ffp, n(i) = A(+) + i(i s 01 ) 

Q, and / = 0,1,2,... the LL numbers. We assume 
sgn(eB) = +. In (0 wc used the E p Ritus' transforma- 
tion E p (x) = 2~^~^±i E pa (x)A(a) (originally developed 
for fermions in |15| and later extended to vector fields in 
[3]), where A(ct) = I±i£^V ; a = ±1, arc the spin 

projectors, and E pa (x) = N n e i( ^ ox0+P2x2+P3X ^D n (p), 
with D n (p) denoting the parabolic cylinder functions 
with argument p = \/1\ell \(x\ —pz/eH), normalization 
factor N n = (AweH) 1 / 4 / Vnl, and positive integer index 
n = n(l,cr) = I + ^2 • can be use d to obtain 

the fermion propagator in momentum space as a function 
of the dynamical mass and the magnetic energy T l 
associated to the AMM of each LL [Io| . 



GW) = (27r) 4 6^(p-p>)Tl(l)G L (p), 



(8) 



with 

&(p) 



0-,CT=±l 



N l ((rT,aV n )-iVj_(A+-Al) 
D l (ooT) 



A(cr)An 



(9) 



and 



A? 




A° 



1 



(1 



2 

-1/ A rl 



■ ivy 



N (aT,aV\\) = aT - M - aV^ , 
D l (aWT) = (M 1 ) 2 - (V[| - oaT 1 ) 2 + (V[) 2 , 
Vf| = (l-Zf | )| P|| | ) Vl = (l-Z l x )\p x \ (10) 



Here, p|| = (p ,0,0,p 3 ) and p ± = (0, 0, v2aBl, 0) denote 
the longitudinal and transverse momenta, respectively, 



3 



while Z'| and Z 1 ^ arc the wave-function renormalization 
coefficients. Taking into account that the MCxSB is an 
infrared phenomenon where the magnetic field becomes 
the leading parameter, we have that the main contri- 
bution to the polarization operator comes from the low- 
energy region where only fermions in the LLL contribute. 
Thus, transforming ([5]) to momentum space and keeping 
the leading contribution I — 0, we obtain 



(<?) 



-2ia\eH\e ^ 



d 2 p 



x Tr 



(2n) 2 

7 llA(+)G (p) 7 llA(+)G°(F^) 



■(H) 



with 



G\p) 



A(+)A+ 



A(+)A^ 



IPllI 



I 



ARA+ 



IP,, | + (M° + T°) 
A(-)A7 



|p„| -(M°-T°) |p||| + (M°-T0) 



(12) 



Note that only the longitudinal components of n ;i „ sur- 
vive in (jlip . Integrating in momenta, using Feynman 
paramctrization, and dimensional regularization yields 



n" (?) 



5 (2),(2) 
: K2" 



(13) 



where 



^2 




(14) 



with £° = M° + T°, the LLL rest energy. Notice that 
the polarization operator ([T3"|) is transverse, q^n^j, = 
n^q" = 0, ensuring the gauge invariance of the LLL 
approximation. In the static limit, go = 0, and infrared 
region |lf | — > 0, the coefficient K2 behaves as 



K2{qo = 0,g-3 -> 0) 



-a|eif| 
3tt(£; ) 2 



(15) 



This result indicates that the inclusion of the AMM term 
does not contribute to produce Debye screening in the in- 
frared region. That is, at distances r > 1/E°, a charge 
within this medium interacts through a normal Coulomb 
potential. Comparing (fT5j) with (j4]) we see that the in- 
duced rest energy E° plays the same role in the broken 
phase of massless QED as the electron mass in massive 
QED. The difference is, however, that E° is not a fixed 
parameter, but it has to be found as the solution of the 
Schwinger-Dyson equation for the electron self-energy. 

To find E° we follow the results of Ref. [Icj . where 
E was found as the solution of the SD equation for the 
electron sclf-cncrgy in the ladder approximation 

E(ac, x') = ie 2 ^G{x, x')j"D^{x - x'), (16) 



where ~E(x,x') is the fcrmion sclf-cncrgy operator, 
D flv {x — x') is the bare photon propagator in the Feyn- 
man's gauge 0, and G(x,x') is the full fermion propa- 
gator {7J. In the LLL, the SD equation reduces to [Hj 



1 = e 2 {AeH) 



d 4 q e-«i 



1 



(2tt) 4 q 2 {E°) 2 +q 2 , 



(17) 



where we introduced the normalized-momentum notation 
q 2 = q 2 /2\eH\. The solution of {17J is given by 



E° ~ J2\eH\exp-J- 
V a 



(18) 



The electric susceptibility can be found now similarly 
to the QED case done previously. From p^|) - ([T5|) and 
(TTBl we have that the medium behaves as a linear, ho- 
mogeneous, and anisotropic dielectric, with electric sus- 
ceptibility in the j^-direction 



XMCxSB 



a\eH\ a 
3tt(E ) 2 = 6^ 



exp 



(19) 



Eq. tfrQj) shows that the susceptibility depends nonper- 
turbativcly on the fine-structure constant and its value 
is independent of the applied magnetic field. Notice the 
marked difference with the QCD situation at strong mag- 
netic field (\eH\ > Aqcd), since in QCD the MCxSB 
leads to a c/iromo-susceptibility that remains a function 
of the magnetic field through the running of the strong 
coupling a s [l?} ■ The colossal susceptibility (fT9|) charac- 
terizes the electric response of the system to an electric 
field parallel to the magnetic one. At zero temperature, 
no critical magnetic field strength is required to catalyze 
the chiral symmetry breaking. On the other hand, even 
though the magnitude of the electric polarization does 
not depend on the magnetic field, for a fixed weak electric 
field, it must satisfies the condition P < ^S)-. That is, in- 
creasing H increases E°, and a larger electric field probe, 
constrained by the condition \eE\ < (E ) 2 < \eH\, is al- 
lowed. We should underline that from a physical point 
of view, it is natural to expect a large anisotropic electric 
susceptibility in the phenomenon of MCxSB in QED, as 
the ground state of the system is characterized by pairs 
forming tiny electric dipoles that can be polarized by the 
external electric field. The role of the magnetic field here 
is to induce the pairs, while the role of the electric field 
is to polarize them. 

The dramatic increase of the electric susceptibility pro- 
duced by the magnetically catalyzed chiral pairs can be 
the best candidate to probe whether the MC\SB mech- 
anism is taking place or not. In a system of massless 
fermions in a magnetic field, if the system exhibits a siz- 
able electric polarization under the application of a weak 
electric field probe along the direction of the magnetic 
field, it will be a plausible evidence of the MC\SB phe- 
nomenon. 

An important implication of this result is that the 
chirally broken phase exhibits strong paraclcctricity, a 
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property found in certain condensed matter systems like 
quantum paraelectric (QP) materials (l8j and transition 
metal oxides (TMO) [19|. In those materials, unaligned 
electric dipoles are aligned in an external electric field, 
producing a high electric susceptibility, often exceeding 
fO 4 . In QP materials the large electric susceptibility is 
temperature-independent below certain critical temper- 
ature, a property attributed to a quantum phase transi- 
tion [lj| ■ An interesting question to explore in the future 
is whether the strong susceptibility found here within 
a (3+l)-dimensional theory is also present in quasipla- 
nar condensed matter systems as bilayer graphene. It is 
known, that the band structure of bilayer graphene can 
be controlled by an applied electric field perpendicular to 
the layers' plane. The electric field creates an electronic 



gap between the valence and conduction bands with en- 
ergy values that varies from zero to midinfrared [20I ] , de- 
pending on the field strength. Under a very weak electric 
field the gap is practically zero and the spectrum is Dirac- 
like. Even though this is a very peculiar 3+1-D system, 
only formed by two layers, one could attempt to model 
it with a 3+1-D theory of interactive massless fermions. 
Because of the universality of the MCxSB, we expect 
that the application of a strong magnetic field parallel to 
the weak electric one will trigger the generation of a dy- 
namical energy gap. Under these conditions, one would 
expect that detecting a very large electric susceptibility 
in the direction of the applied fields would signal the re- 
alization of the MCxSB mechanism. 
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